Higher-order defect mode laser 
in an optically thick photonic crystal slab 
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The use of an optically thick slab may provide versatile solutions for the realization of a current injection type 
laser using photonic crystals. Here, we show that a transversely higher-order defect mode can be designed to 
be confined by a photonic band gap in such a thick slab. Using simulations, we show that a high-Q of > 10 is 
possible from a finely tuned second-order hexapolc mode. Experimentally, we achieve optically pumped pulsed 
lasing at 1347 nm from the second-order hexapole mode with a peak threshold pump power of 88 /iW. © 
2012 Optical Society of America 
OCIS codes: 230.5298, 250.5960. 
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Two-dimensional (2-D) photonic crystal (PhC) slab 
structures have, so far, been in the form of a thin di- 
electric slab, whose thickness T is often chosen to be 
^200 nm for an operational wavelength of ~1.3 /im. 
This thickness consideration is to maximize the size of 
the photonic bandgap (PBG) in the in-plane direction 
(x-y plane) [1], which has unfortunately placed a severe 
constraint on the design of a current-injection type laser. 
Pulsed lasing operation has been demonstrated using a 
vertically- varying p-i-n doping structure within the thin 
PhC slab, for which a sub-micron size dielectric post 
placed directly underneath the laser cavity serves as 
a current path [2]. Recent efforts have moved towards 
a laterally-vaxy'uxg p-i-n structure and a few successful 
results were already reported by groups in both Stan- 
ford [3] and NTT [4]. However, there are still favorable 
reasons for using a vertically-varying doped structure, 
because such a design allows a monolithic growth of all 
of the epitaxial layers that are almost free of crystal de- 
fects. 

Recently, we have shown that even a very thick slab 
can support sufficiently high-Q (few thousands) cavity 
modes for lasing. [5] In our previous result, however, the 
dipolc mode formed in a triangular lattice air-hole PhC 
slab was emitting more photons into the in-plane direc- 
tions rather than into the vertical direction (z) for effi- 
cient photon emission and collection. Moreover, Q could 
not exceed 3,000 with T = 606 nm. It would seem, at 
first, that we have no other options for further improve- 
ment in Q, since the poor horizontal confinement appears 
inevitable due to the absence of a PBG. It is our purpose 
in this Letter to rebut this first intuition and show that 
the thick slab can be used to achieve an efficient vertical 
emitter with a surprisingly high Q of over 10 5 . 

To start, we perform numerical simulations both us- 
ing the plane- wave-expansion method (PWE) [6] and the 
finite-difference time-domain method (FDTD) to inves- 
tigate how a PBG evolves as we change the air-hole ra- 



dius (r) and the slab thickness (T) [Fig. 1(a)]. Note that 
r and T are represented in the unit of the lattice con- 
stant (a). In the case of a triangular lattice air-hole PhC, 
{r = 0.40a, T = 0.6a} gives the widest gap centered at 
uj c « 0.38, which agrees with earlier work by Johnson, 
et al. [1]. Also note that there exists a broad region of 
{r, T} that gives a wide gap-to-midgap ratio [1] Alu > 
30%. This is why r and T are often chosen to be ~ 0.35a 
and 0.5a, respectively. We also find that a tiny PBG 
(usually Aw ~ 1%) exists up to T = 1.25a. 

The dipole mode discussed in our previous work [5] 
(T = 1.86a) is marked as 'Id' in the gap map. Now, 
we pose the question of whether we can design a certain 
resonant mode emitting at ~1.3 fim that is confined by 
a PBG in a slab with T = 606 nm. From the gap map 
diagram, the only possibility appears to be increasing 
a in order to bring down T(a) below 1.25a. However, 
keeping the same 'Id' mode, larger a usually results in 
the longer A, because u> — a/X is rather fixed by the in- 
planc modal structure of a resonant mode [7]. Therefore, 
we should look instead into other resonant modes that 
do not resemble the dipole mode. 

It is well known that even a single defect resonator 
supports multiple resonances such as the quadrupole, 
the hexapole, and the monopole modes [8]. These higher- 
order modes are pulled down from the conduction band- 
edge of the photonic band structure [7] . Further tuning 
the defect region can get more higher-order modes pulled 
down into the gap. One possible route from the (first- 
order) dipole mode ('Id') to the second-order hexapole 
mode ('2h') is drawn by an arrow in the gap map di- 
agram. The '2h' is designed to be resonant at a wave- 
length close to that of 'Id' [9] even though it has quite 
a large a of 500 nm (thus, T « 1.21a). As a quantita- 
tive measure showing how well the PhC layers work as 
a mirror, we calculate the vertical extraction efficiency 

?7 vcr t defined by ?7 vcr t = (l/Qvert)/(l/Qhorz + VQvert) = 
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Fig. 1. (a) A 2-D map of a PBG for a triagular lattice air-hole (radius = r) PhC in a dielectric slab (n s i a b = 3. 4) with 
a thickness of T. The 2-D color scale map represents the size of the PBG in terms of the gap-midgap ratio defined 
by Aw = Aw/w c , where w c is the center frequency of a PBG. The contour lines of w c are overlaid on the 2-D map. 
Note that throughout the Letter, all frequencies are normalized by 2irc/a, hence w = a/A (dimensionless) . (b) The 
(first-order) dipole mode [Q=2,600 and V = 0.82(A/n s i a b) 3 ] oscillating at A = 1341 nm with a = 325 nm and (c) the 
second-order hexapole mode [Q=15,200 and V = 2.23(A/n s i a b) 3 ] oscillating at A = 1365 nm with a = 500 nm. Both 
modes arc formed in a slab with T = 606 nm. 



(l/Qvcrt)/(l/Qtot) [5]. We find that rj VCTt of '2h' shown 
in Fig. 1(c) is 0.954 (Q h orz = 3.3 x 10 5 ) with the same 
number of air-hole barriers shown in Fig. 2(b). We be- 
lieve this ?7vert (or Qhorz) has not yet been saturated due 
to the small gap size, expecting further improvement by 
increasing the number of barriers. Probably, in apply- 
ing the idea of a higher-order resonant mode, T of 606 
nm would be the upper limit for A ~T300 nm, as '2h' 
can only be made barely located at the top-right corner 
of the gap map diagram. We would like to note that the 
same strategy can be applied more effectively to the case 
of an intermediate thickness range of 400 nm < T < 600 



(a) 



4< 

,' t — r 

Ri R2 R3 R 4 R5 R>( 




R, 


= 0.26x« 




= 1.07xrt 


K 2 


= 1 ,04x« 



R 2 = R - 0.08xfl 
R 3 = R - 0.04xrt 
R 4 =R 

R 5 =(R+R b0 )/2 

R>6 = Rba 



JOC 

000000 

)000000 
QQQPOOO 

3O0QOOOO 

QOQQpOOO 
OOOOOOOOO 
^OQQOOOO 
^iPQOOOOOO 
#^)QDOOOOO 

QQQOOOOO 
^OOOOOOO"" 
OOOOOOO 
3000000 
OOOOOO 
DOOOOO 
OOOOO 



OOO 
)OO max , 

00 

)0 IEI 2 

o 



Fig. 2. (a) A schematic diagram shows how we finely 
tune air-hole sizes and locations to optimize Q. (b) An 
electric-field intensity distribution (|E| 2 ) of the highest- 
Q mode (case II in Table 1). 



nm. Imagine a first-order resonant mode ('lx') oscillat- 
ing at uj w 0.26 within a slab with T = 1.1a. At this 
region, Aw is only about 5%. We can bring it down deep 
into the band gap by utilizing its second-order resonant 
mode ('2x'). If '2x' oscillates at ui m 0.33, then, without 
altering A, '2x' can be formed in a slab with T pa 0.87a, 
at which Aw is as large as 20%. 

Surrounding the '2h' with the large air-holes of R = 
0.46a would give better spectral matching between the 
'2h' resonance and the center of the tiny bandgap. How- 
ever, such a large air-hole radius is not advantageous for 
the device's mechanical robustness. Therefore, we pro- 
ceed to study if the background air-hole radii (Rbg) can 
be substantially reduced without sacrificing Q too much. 
Several representative cases of fine-tuned air-holes are 
shown in Fig. 2 and Table 1. u) and Q are most sensi- 
tively dependent on the parameters near the center of 
the resonator; R\, K\, and K2- These parameters have 
been determined in a manner to optimize Q. The out- 
skirt region from R4 is intended as a mirror. Air-hole 



Table 1. Examples of the second-order hexapole 
mode in a T = 606 nm slab. 



case 


R(a) 


-Rbg (a) 


Qtot 


Qvert 




I 


0.45 


0.45 


55,400 


58,500 


0.947 


II 


0.45 


0.38 


105,100 


146,200 


0.719 


III 


0.44 


0.38 


50,400 


63,900 


0.789 


IV 


0.43 


0.38 


27,900 


34,400 


0.811 


V 


0.42 


0.38 


17,800 


21,800 


0.813 


VI 


0.41 


0.38 


12,400 


15,300 


0.807 
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Fig. 3. (a) Scanning electron microscope (SEM) image 
taken at a tilt of about 10°. Note that T = 606 nm. (b) 
FDTD simulated mode profile. The actual SEM image 
was used to directly input air-hole shapes and sizes, (c) 
Light-in versus light-out (L-L) curve. Insets show near- 
infrared camera images taken before and after the lasing. 
(d) Lasing spectrum measured at a peak pump power of 
200 fxW. 



radii before and after R4 are designed to vary gradually 
to minimize unintentional scattering losses at the crystal 
dislocations. Since R\, K\, and if 2 are fixed, all the res- 
onant wavelengths tend to stay near 1323 nm. a = 450 
nm for all those cases, thus T = 1.35a and there exists 
no PEG. 

Contrary to the initial expectation, Q can be made 
higher even in the absence of a rigorous PEG [10]. In 
Case II, we find that Qvert can be greatly improved by 
more than a factor of 10, thereby Qtot can reach over 10 5 . 
It is interesting to observe that, comparing I and II, air- 
holes located far from the mode's energy (i?>6) can affect 
Qvcrt- It should also be noted that just one layer of R4 = 
0.45a effectively blocks the horizontal photon leakage. As 
we progressively reduce R, both Q V ert and Qhorz decrease 
somewhat. At the final stage of the tuning (VI), all air- 
hole sizes become reasonable for experimental realization 
and Q remains well above what is required for lasing. 

In experiment, we intend to fabricate structurally 
more robust design similar to VI rather than the Q- 
optimized design of II. We use the same InGaAsP wafer 
containing 7 InGaAsP quantum wells emitting near 1325 
nm used in our previous work [5] . To define high-aspect 
ratio air-holes, we use chemically-assisted ion-beam etch- 
ing with Ar and CI2 [5] . The fabricated devices are opti- 
cally pumped at room-temperature with a 830 nm laser 
diode driven by a pulse generator at 1 MHz with a duty 
cycle of 2.5 %. A 100 x objective lens is used to focus 



the pump laser on the center of the resonator. The L- 
L curve clearly shows a threshold, estimated to be 88 
fjiW in terms of peak pump power, where we have as- 
sumed about 20% of actual incident pump power is ab- 
sorbed by the slab. We verify single mode lasing oper- 
ation over a wide spectral range (1300 nm~T400 nm) 
with a side-mode suppression ratio of ^30 dB. To con- 
firm if the measured laser peak truly originates from the 
'2h', we perform FDTD simulation using a contour in- 
put for actual fabricated air-holes from the SEM image. 
The FDTD expects that the designed '2h' mode should 
locate at a wavelength of 1340 nm, which agrees very 
well with the experimental result. 

In summary, we show that a PhC slab with optically 
thick T — 606 nm can be used to construct a PBG- 
confincd resonant mode oscillating at a wavelength of 
^1300 nm. We also show that a surprisingly high Q of 
over 10 5 can be obtained even in the absence of a rigorous 
PBG, and that the use of the higher-order resonant mode 
can be quite advantageous for making an efficient PhC 
laser with an optically thick slab. 
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